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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS — 
TECHNICAL MEMORANDUM NO. 1349 


ON A CLASS OF EXACT SOLUTIONS OF THE EQUATIONS OF MOTION 
OF A VISCOUS FLUID* 


By V. I. Yatseyev 


The general solution is obtained herein of the equations of motion 
of a viscous fluid in which the velocity field is inversely proportional 
to the distance from a certain point. Some particular cases of such 
motion are investigated. 


1. The motion of a viscous fluid with velocity field and pressure 
in spherical coordinates can be given by the following expressions: 


„hl = = Ye =9 5= els) (1) 


A particular solution of the equations of Navier-Stokes for this case 
was obtained by Landau (reference 1). In the present paper a general 
solution is given of the equations of Navier-Stokes for the motion of 
the class under consideration. 


Substituting expressions (1) in the equations of Navier-Stokes and in 

the equation of continuity yields the following system: 
Fo + fÊ - fF' + 2g -V E + F' cot 0 - 2f' - 2F - 2f cot o] =0 
(2) 


ES cot 8 + 2F' - f (1 + cot? )| = o (3) 


F+f' +f cot 6=0 (4) 


Determining F from equation (4) and substituting in equations (2) and 
(3) give 
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£12 + ff" + 3ff' cot 0 + 2g -V Es + 2P" cot 8 - f' (2 + cot’ 9) + 


f cot 9 (1 + cot2 o)| =0 (5) 
ff’ +8" +Y E + f' cot @- f (14 cot? 9) =0 (6) 


Differentiating expression (6) 


P'E + ff" + g" +d E + f" cot 0 - 2f' (1 + cot? 0) + 


2f cot 8 (1 + cot® o) | 20 (7) 


Eliminating the nonlinear terms fre, ff", and ff' from equation (5) 
with the aid of equations (6) and (7) yields a linear equation in the 
function g + avf': 


(g + 20f')" + 3 cot 9 (g + 20£')' - 2 (g + 20f') =0 (8) 


the general solution of which is in the form 


g+ avr = 2p? bcos 0 ~ a (9) 
sin © 


where 2V2a and 202b are constants of integration. 


Integrating equation (6) 


f2 + 2g + 2V (f' + f cot 0) = - 2p%c (10) 
where 20%c is the constant of integration. 


The function g(@) is eliminated from equations (9) and (10) to 
give an equation of the Riccati type for the function f:l 


f! =g fË +f cot o + zo| (11) 


lafter sending the manuscript to press the author obtained from 
L. D. Landau a communication on the work of N. Slezkin (reference 2) in 
which he arrived at the same equation by a different method. 
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The substitution 


f = - zux'(0)/X(e) (12) 


reduces equation (11) to the linear equation: 


x" -x' cot 0 + b con f san s x=0 (13) 
sin” 0 


which by the substitution 


z = cos? (9/2) (14) 


is transformed into an equation of the Fuchsian type: 


dx a+rb-2(b+c)z + Zeze 


=0 15) 
dz? 4z° (z - 1) x ( 


The usual computations (reference 3), which are omitted herein, give the 
general solution of eguation (15) as: 


l+atß- 
x(0) = (cos 2) (se z) is cıF (or cos? 2)* 


coF (a + 1-Yv,B+1-y,2- y, cos“ 2) (16) 


where the parameters of the hypergeometric function a,B,y (which can 
also have complex values) are connected with the constants of integra- 
tion a,b,c by the formulas: 


2 
a=y?-(1+0+8) y+ fo + 6)" 


ya 
2 


A CITA (17) 


O 
il 
DH 


c= la- B)2-1 
2 


4 NACA TM 1349 


Formulas (4), (9), (16), and (17) give the general solution, depend- 
ing on the four constants a, b, c, and A= co/cy of the Navier-Stokes 


equations for the class of motion of a viscous fluid under consideration. 
The constants of integration a, b, and c are expressed in terms of the 
corresponding tensor components of the density of the momentum transfer: 


Ilik = Põik + PViVk - OV —— + — 18 
i ik ivk Sub: zi 
Carrying out the computations 
2 
1 _ 2V“p fb cos 0 -a 
PP rê sin? 0 
2vĉp a - b cos 9 c 
le were (19) 
r sin” 0 
a 2vĉp c cos 9 -b 
Mae RE 
r sin 0 
The streamlines are determined by the equation: 
dr/v, = rad0/vg (20) 
the integration of which gives 
const/r = f sin 9 (21) 


2. Attention is now given to two particular examples for which the 
equation of Fuchs degenerates. 


(a) Equation (15) has only one regular singular point, z = œ. In 
this case 


a=b=c=0 (22) 


and therefore by equations (19) 


Topp = TIgg = Tirgo =0 (23) 
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The particular solution of equation (15) 


KOO) =2z -1-A (24) 


leads by formulas (4), (9), and (12) to the solution found by Landau: 


2 f 
A” - 1 2v sin 9 
= A o = —————— 
ne i a - cos 8)? A me) en 
g(@) = 4v2 l-A cos 0 (25) 


(cos 9 - A)? 
This solution is analogous to the problem of a stream flowing out of the 
end of a thin pipe into a region filled with the same fluid. It is the 
only regular solution for all values of the angle 9. 


(b) Equation (15) has only two regular points z =0 and z =o. 
In this case it follows from equation (15) that 


a=b=c #0 (26) 
and equation (11) becomes Euler's equation 
222 (a /az2) - ax =0 (27) 


the general solutions of which are 


x(0) = ex/2 cosh (nx + A) fora> - 1/2 
x(0) = ex/2 cos (nx + A) fora< - 1/2 (28) 


x(0) ex/2 (1 + Ax) for a= - 1/2 


Correspondingly, the following equations are obtained for the function 
£(9): 


H 
li 


2v => (m tanh (nx + A) + ya for a > - 1/2 (29) 


sin 8 1 
= u nF <- 
f = 20 I 5 £ - n tan (nx + wh for a 1/2 (30) 


f= 2v Bids (te + ve) for a= - 1/2 (31) 
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where -x = ln (1 + cos 8), n= 5 YA + Za |. (For a =0, n = 1/2 in 
equation (29) the solution of Landau is again obtained. ) 


For the solution of equation (29) by formula (4) 


l - cos 8 nê 


l + cos 0 Sosh” (nx + A) 


(32) 


F(9) = - 2v {n tanh (nx + A) + 1/2) + 2V 


while g(@) is determined by formula (9). 


The equation of the streamlines is in the form 


1 
const/r = (1 - cos 0) {m tanh (nx + A) + WA (33) 


where the values of the constants n and A are determined from the 
conditions 


r(3)- 20 (n tanh A + 1/2) 


: (34) 


n K n 
r(5)+* (5)- E A SAA 
2 2 cosh? A 


The obtained solution corresponds to the problem of the stream flowing 
from the half line 9 = x into a region filled with the same fluid. 


For solution (30), the parametric equation for the streamlines is 


in the form 


const/r = (2 - eX) EZ - n tan (nx + a) | 


9 = arc cos (eë - 1) (35) 


The function (35) for 6x isa strongly oscillating one. It can 
therefore be concluded that solution (30) has no physical sense. 


The author acknowledges the suggestions and help received from 


Professor Y. B. Rumer. 


Translated by S. Reiss 
National Advisory Committee for Aeronautics 
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